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Abstract 



We prove the following results, (i) One-dimensional Bose gases which interact via unsealed inte- 
grable pair interactions and are confined in an external potential increasing faster than quadrat- 
ically undergo a complete generalized Bose-Einstein condensation (BEC) at any temperature, in 
the sense that a macroscopic number of particles are distributed on a o(N) number of one-particle 
states, (ii) In a one dimensional harmonic trap the replacement of the oscillator frequency u by 
cj In N/N gives rise to a phase transition at a = hu/3 — 1 in the noninteracting gas. For a < 1 the 
limit distribution of no/N a is exponential and {no)/N a — > 1. For a > 1 there is BEC with a con- 
densate density (no) /AT — » 1 — a -1 . For a > 1, (In N/N) (no — (no)) is asymptotically distributed 
following Gumbel's law. For any a > the free energy is — {-k 2 /G(3a)N / In N + o(N/]nN), with 
no singularity at a — 1. (iii) In Model (ii) both above and below the critical temperature the gas 
undergoes a complete generalized BEC, thus providing a coexistence of ordinary and generalized 
condensates below the critical point, (iv) Adding an interaction (Un) = o(N\nN) to Model (ii) 
we prove that a complete generalized BEC occurs for any (3 > 0. 

PACS numbers: 0530J, 0550, 7510J 

KEY WORDS: trapped Bose gas; one dimension; generalised Bose-Einstein condensation; in- 
teractions 

1 Introduction 

The idea of Bose-Einstein condensation (BEC) in the ground state of a one-dimensional interacting 
Bose gas first appeared in Girardeau's 1960 paper about the <5-gas in the impenetrable limit [1]. Based 
on an approximate calculation, Girardeau suggested that there is no BEC in the ground state but, 
instead, there is a complete generalized BEC (GBEC) in the sense that 



Here N is the number of particles, (n&) is the ground state expectation value of the occupation number 
of the one-particle state ~ exp(ikx) and p is the density. Later Schultz [2] disproved this conjecture by 
showing that the double limit in (1) yields zero. A subsequent work by Lenard [3] on the momentum 
distribution confirmed this conclusion. Somewhat later the absence of BEC [4] (see [5] for a rigorous 
proof) and that of GBEC [6] were shown in one- and two-dimensional homogenous systems at positive 
temperatures T. Bogoliubov's inequality [7] which these results are based on becomes trivial at T = 0, 
and BEC or GBEC in the ground state of the soft- (5-gas has been a most interesting open problem 
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ever since the Bethe-Ansatz solution by Lieb and Liniger was published [8]. Different approximate 
methods, as an effective long-range theory [9], expansion to first order in 1/c where c is the prefactor 
of S [10] and conformal field theoric approach [11], predict an algebraic decay of the ground state 
expectation value of the boson field operator a*(x)a(0), incompatible with an off-diagonal long-range 
order. Also, some arguments [12, 13] excluding the spontaneous breakdown of a continuous symmetry 
in the ground state of certain one-dimensional systems probably apply to the 5-gas. On the other 
hand, no conjecture seems to exist concerning the possibility of GBEC. 

BEC can occur in one dimension if the gas is confined in an external potential. In [14] it was 
shown that there is BEC at any finite temperature if the bosons are in a fixed (not TV-dependent) 
external potential and interact via a suitably scaled pair interaction, 

u N (x) = b N u(a N x) . (2) 

Here u is an integrable positive function (a soft 5 in one dimension is allowed) and 6 at and a at have 
to be chosen so that 

J u N dx < C/N (3) 

for some finite C. This result is valid in any dimension d > 1 and imposes 

b N < Ca d N /N. (4) 

Interactions satisfying this condition cover in any dimension a whole range from mean- field types (ajv 
and, thus, bN tending to zero as N goes to infinity) to sharply concentrated ones (un and 6 at going 
to infinity). Even the mean-field type interactions are nontrivial if the confining potential is locally 
bounded, thus allowing particles to decrease their interaction by increasing their separation, or if u(x) 
tends to infinity as x goes to zero. 

For the noninteracting gas Bose condensation in a fixed trap is a pure ground state phenomenon. 
Let H° = — t^A + V where V is chosen in such a way that tre _/3H < oo for any (3 > 0. (By 
Symanzik's version [15] of the Golden-Thompson-Symanzik inequality 

tr e -^° <\- d Je-^dr (5) 

where A/3 = hy / 2Tr/3/m, this holds if V increases faster than logarithmically.) As it was shown in 
[14], at any temperature the N — > 00 limit of the thermal equilibrium state of the corresponding 
noninteracting gas remains essentially a finite perturbation of the ground state. For this reason, at 
any finite temperature there is an asymptotically complete condensation into <po, the ground state of 
H°. Because the ground state is common for the Bose and Boltzmann gases, there is BEC, although 
incomplete, at any temperature even in the trapped ideal quantum Boltzmann gas! The reduced one- 
particle density matrix of the latter is Ne~P H / tier^ H , whose largest eigenvalue JVe"^ / tre _/3ff 
gives the mean number of particles in ip (e is the lowest eigenvalue of H ). 

Adding a scaled pair interaction (2) with property (4) to the noninteracting gas while keeping 
the trap fixed has a nontrivial effect, which can be expected by comparing energies: The total free 
energy of the ideal gas is of the order of 1 (apart from the trivial ground state energy Nsq which 
could be chosen to be zero), while the total interaction energy under condition (3) is of the order of 
N. BEC survives such a perturbation at all temperatures. Working directly at zero temperature, in 
three dimensions for bN = N 2 and cxn = N (which is Gross-Pitaevskii scaling, satisfying (3)) Lieb 
and Seiringer were able to prove a complete BEC into the minimizer of the Gross-Pitaevskii energy 
functional [16]. This is very different from ip , also showing the nontrivial effect of the interaction. 

The strength of an integrable interaction can be measured either by its integral or by its scattering 
length. Whether or not a pair interaction un whose integral is of order 1/N counts to be weak from the 
other point of view, depends on the space dimension. If wat(x) = a N u(aNx), the scattering length of 
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un is 1/ajv times the scattering length of u. In three dimensions (3) imposes oin ~ N, thus both the 
integral and the scattering length are of order 1/N. In two dimensions J a N u(aN^) dx = J u(x) dx, 
so un has to be chosen in a different form in order to comply with (3). In one dimension, with the 
necessary choice ctN ~ 1/N, we get a scattering length ~ N, therefore the two criteria of strength 
contradict each other. This last example is a mean-field type interaction while e.g. a at = N 2 and 
&at = N represent the opposite limit in one dimension. 

Further insight is obtained if scaling of the interaction is transformed into scaling of the potential 
and the temperature: 

N h 2 

i=l \<i<j<N 

[ N h 2 ) 

= a 2 N lY,{-^ A y> + a N 2 v(aN 1 y l )}+ E u (y*-yi)\ ( 6 ) 

[i=l l<i<j<JV J 

where y» = chatx^ Because a 2 N multiplies the inverse temperature (3, we obtain a joint limit during 
which N — > oo and in three dimensions the temperature goes to zero and the trap opens in such a way 
that the particle density tends to zero, while in one dimension the temperature goes to infinity and 
the trap closes so that the density diverges. The three dimensional example suggests that whatever 
high the unsealed temperature, Gross-Pitaevskii scaling ajv = N may reduce the thermal equilibrium 
state to the ground state. It is at least true that BEC remains complete at positive temperatures (see 
[17] for the Gross-Pitaevskii limit and the discussion of Section 2 for the general case (2)). 

The analysis of the behaviour of a many-body system in a scaled external field or in the presence 
of scaled interactions comes about naturally in the study of trapped Bose gases. Taking the limit 
TV — > oo and scaling certain quantities during this limit is an abstraction to obtain a mathematically 
clear-cut answer, not different in its philosophy from the thermodynamic limit in homogenous systems. 
The most frequently used scaling limit is related to the Gross-Pitaevskii energy functional [18, 19] 
which describes the ground-state properties of dilute Bose gases. As mentioned earlier, in a three- 
dimensional trap this scaling is characterized by fixing NIn, where In is the scattering length of the 
pair interaction, see also [16]. Considering the scattering length as an effective diameter of the particle, 
space filling is seen to decrease as 1/N 2 in this limit. Condition (2)-(4) is a generalization of the 
Gross-Pitaevskii scaling of the interaction. As a different example, scaling of the external confining 
potential appears in the theoretical discussion of BEC in elongated traps, realized experimentally 
more recently [20, 21]. Various aspects of a one-dimensional behaviour can become manifest because 
a tight transverse trapping makes transversal degrees of freedom freeze out. Theoretically this can be 
achieved by applying different scalings of the confining potential in the transverse and longitudinal 
directions. These studies reveal a rich structure of low-dimensional regimes [22, 23, 24]. Note also 
that a theoretical work on one-dimensional bosons in a scaled harmonic trap was published well before 
the first experimental realization [25] . 

The present paper is motivated by two naturally arising questions. What happens with BEC in 
dense or strongly interacting trapped gases? How does interaction modify BEC occurring with a phase 
transition in the ideal trapped gas? To study these questions, first (Section 2) for a gas in a fixed trap 
we relax the condition (3) at the price of not being able to prove BEC, only GBEC. We obtain the best 
result for one dimensional anharmonic and other supcrharmonic traps, i.e. potentials with a faster 
than quadratic increase. We can prove a complete GBEC at any finite temperature in the presence 
of unsealed pair interactions, that is, the total interaction energy increasing as TV 2 . This result holds 
also for attractive interactions, showing that the collapse of an attractive Bose gas can be considered 
as a GBEC. A box is a superharmonic trap and thus, in principle, our findings may have implications 
for the homogenous 5-gas. However, we can prove GBEC only if c/p — > as N — > oo. In other words, 
we just fail to prove it for the interacting gas. A reason of this failure may be that probably there 
is no Bose condensation, ordinary or generalized, in the ground state of the homogenous soft-5-gas. 
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Although our method is too weak to prove BEC in homogenous interacting gases, some nontrivial 
bounds can be obtained and are presented in Section 2 (equations (29) and (30)). 

In the case of a gas in a fixed trap our proof either yields BEC/GBEC for any finite temperature 
or does not yield it at all. In the second part (Section 3) of the paper we are interested in the stability 
against interactions of BEC which occurs through a phase transition in the noninteracting gas. Because 
T c = oo for fixed traps, it is natural to look for a finite T c by opening the trap together with N — > oo. 
During this limit the separation of energy levels and, in particular, the gap above the ground state tend 
to zero, a situation similar to that occurring in the homogenous gas. As a consequence, the difficulties 
to include interactions start to resemble those appearing in homogenous systems. We consider only the 
harmonic potential, V(x) = mu> 2 x 2 /2. In this case the relevant dimcnsionless parameter is a = Hlj/3. 
In three dimensions it was found [18] that the iV-dependent critical temperature of the noninteracting 
gas is given by kBT c (N) = fuvlN/^fi)] 1 / 3 which is equivalent to say that the replacement of a by 
aN ' gives rise to BEC at a c = C(3) 1/3 « 1.063. Although the one-particle density of states 
is different, this scaling qualitatively reproduces the properties of the homogenous gas. Since the 
interactions for which we can prove the survival of condensation are too weak in three dimensions 
((^(iV 1 / 3 ), negligible on the scale N of the free energy of the noninteracting gas), in Section 3.1 wc 
consider the analogous problem in one dimension. Earlier Ketterle and van Druten [25] found that in 
the noninteracting gas the iV-dependent critical temperature was IcbT c (N) = fkuN/ \n(2N).We start 
with a detailed study of the noninteracting system. We replace u> by lo"/n or a by ajw where 7jv 
tends to zero and discuss the different possibilities. Depending on the decay rate of 7jv, T c may be 
infinite or zero or may have a finite positive value. If jn = In N/N, T c is finite positive with a c = 1 
(equivalently, k B T c (N) = hwN/]nN; although wc note that T C (N) is ill-defined and 7^ = \naN/N 
leads to a c = 1 for any fixed a > 0) and with BEC for a > 1 and no BEC for a < 1. On the other 
hand, T c = 00 if j N /(lnN/N) -> 00 and T c = if 7;v/(ln N/N) -> 0. In particular, for 7^ = l/N 
we find T c = and an extensive free energy and can, therefore, identify this scaling limit with that 
of the homogenous system. Since exact computations are possible, for the critical scaling we can 
obtain the asymptotic distribution of n , the number of particles in the ground state, and its mean 
value together with finite-size corrections. We find that for a < 1 the limit distribution of no/N a 
is exponential and (n )/N a — > 1. For a > 1 the condensate density (no)/N — > 1 — a -1 and the 
fluctuation of n is huge, (In N/N)(n — (n )) is asymptotically distributed following Gumbel's law. 
This holds true at the critical point a = 1 as well where (n ) = N In In N/ In N + 0(N/ In N). Not 
surprisingly, the free energy is subcxtensive, = — (n 2 /Q[3a)N/ In N + o(N/lnN) for any a > 0, 
with no singularity at a = 1. A subtle difference compared with the three dimensional case, revealing 
itself only through the details described above, is that at all temperatures the gas forms a complete 
generalized Bose-Einstein condensate with the participation of at most <~ N/ In N one-particle levels. 
We prove this in Section 3.2. We also show there that for a < 1 a large number of low- lying levels arc 
equally occupied by iV a particles, and for a > 1 particles not condensed into the one-particle ground 
state are in a generalized condensate at its critical point and there is no fragmented condensation, 
(rii) — o(N) for each i > 0. In Section 3.3 we study the effect of interactions. We can only prove that 
at least a complete GBEC occurs at all temperatures provided that {Un)i3h /NlnN goes to zero 
as N tends to infinity. Here Un is the iV-particle interaction energy and (Un)/3h% i s its mean value 
taken in the thermal equilibrium state of the noninteracting gas in the scaled trap. Thus {Un)(3h° n 
can be of the order of N or of slightly higher order, so that F^/{Un)ph° n ~ * 0j showing that the effect 
of the interaction is indeed nontrivial. For scaled pair interactions the analog of (4) is the somewhat 
weaker condition (stronger interaction) = o(ajv IniV/iV) . The mathematical analysis in Section 
3 is conceptually not difficult but it is ramified and lengthy. Readers uninterested in details can go 
directly to the assertions of Theorems 3, 4 and 5. The paper ends with a Summary. 
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2 Generalized Bose condensation in traps 



2.1 Preliminaries 

The Hamiltonian of N nonintcracting bosons is 



< = Et-^ A - + V (*i)} = ®f=i^° • (7) 

i— 1 

We allow the confining potential V and thus the one-particle Hamiltonian H° to depend on N, and 
omit to indicate the resulting TV-dependence of the eigenvalues £o < £i < £2 < • • • and eigenvectors 
ipj. The iV-particle interaction energy and the corresponding Hamiltonian will be denoted by Un 
and Hn = H N + Un, respectively. Un is supposed to be bounded from below. As in [14], for a 
positive integer J and a positive number 5 we define the modified Hamiltonians H°(J,5), H N (J,5) 
and Hn(J, 5) by replacing Sj with Sj + 6 for j < J in the spectral resolution of H°. For an operator 
A, (A)f)H = Tr A e~P H / Tr e~@ H with the trace taken in the iV-particle symmetric subspacc. In 
particular, (rij)pH is the mean occupation number of ipj. 

The basic estimate we use for proving BEC or GBEC in the interacting gas is the following. 

Lemma 2.1 For any [3 > 

J J ^ 

J2( n j)fJH N > J2(n 3 )pH%(.LS) - -g[(UN}f)H° N -infC/jv] 
3=0 j=0 



> 

This bound is a consequence of the identity 



n o)pH%(j,8) - \ [{Un)ph° -inf U N ] . (8) 



.7 o 

2Lj n 3/0H^(J,S) d(f35) ' [) 

the convexity of the logarithm of the trace as a function of (38 and Bogoliubov's convexity inequality 
[26] 

TV p—PH-L 

cf. the first part of the proof of Theorem 4.2 in [14]. 

Lemma 2.1 is useless in the case of hard-core interactions which yield —00 on the right-hand side 
of the inequality (8), but in the case of integrable interactions it leads to nontrivial results. Suppose 
that V is fixed, and 

t/jy(xi, ■ ■ ■ ,Xjy) = y^Ujv(Xt - Xj) (11) 

i<3 

One can show that 

^K^ivWo - (*o,EWI < c(/3)\\un\\i (12) 

where $0 = ¥?o( x i) ••• yo( x Jv)j 1 1 tijv 1 1 1 = / l u iv(x)| dx and c(/3) is independent of iV. The above lemma 
and the bound (12) were used earlier to prove a theorem (Theorem 4.2 of Ref. [14]) that we repeat 
here in a slightly different form. 

Theorem 1 Suppose that V is independent of N . Let Un be a positive pair interaction with ||mjv||i < 
C/N for some constant C. Then 

L(U) = lim ($ , U N $ )/N < 00 , (13) 

JV^oo 
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and for any (3 > and J > 

J 

>tfT ooN v>"""''-- e J+1 -e 



E lim W > 1 - • (14) 

The main example is (2) with (4). Because L(U) is hnitc and ej are independent of N and tend to 
infinity, there exists a Jo independent of N and uniquely defined through 

£.Jo ~ £o < £(£0 < £Jo+i - £o • (15) 

The right-hand side of (14) is positive for J > Jo, implying that at least one of ifo . . . , ip,j is micro- 
scopically occupied. This is just Bose-Einstein condensation, because of the following simple result, 
shown in [14]. 

Lemma 2.2 Let a(N) be any N -particle density matrix and ip any normalized element of the one- 
particle Hilbert space H. Then 

(n[ip]) a = (ip,(7iip) (16) 

where oi(iV) is the one-particle reduced density matrix corresponding to a and n[ip] is the occupation 
number operator in H N , associated to ip, 

n[ip] = \ip){ip\® I ® ■ ■ ■ ® I + ■ ■ ■ + I ® ■ ■ ■ ® I ®\ip){ip\ . (17) 

Thus (<Pj,(?i<Pj) — {n[pj\)i3H N = {nj)f3H N is of the order of TV for at least one of j = 0, . . . , Jo, and 
the maximum eigenvalue of u\ satisfies 

M > » (i-^U, (i 8 , 



Jo + 1 V £./o + l _ £ 

implying BEC. Although in [14] the theorem was proved for a general J as presented above, it was 
stated only for J = Jo, because we were interested in finding the minimum number of eigenstates of 
one of which can be seen to be occupied macroscopically. In fact, Theorem 1 implies a complete 
BEC at all temperatures. To see this, it suffices to take the limit J — > oo in (14), yielding 

1 

3=0 

Only (nj)pH N oc N contribute to the sum. Thus, all but an asymptotically vanishing fraction of 
particles are carried by macroscopically occupied one-particle eigenstates. It is not difficult to show 
that the same is true for the eigenstates of a\. If a\ipj = \ji>j, j = 0, 1 . . ., 

oo ^ oo 

E J im m <»fe]w = E J im x i ( N )/ N = 1 ( 2 °) 

3=0 3=0 

comes from (14) and the following generalization of the variational principle. 

Lemma 2.3 Let A be an upper semibounded self-adjoint operator on a separable Hilbert space H. 
Suppose that A has a pure point spectrum Ai > A 2 > . . .. Then for any positive integer n 

n n 

E A *= su p E(^'^)- ( 21 ) 

i=l 01, . . . ,<f> n € H i=l 

{(j>i,<j>j) = Sij 
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Proof. Let tpi be the orthonormal eigenvectors of A. Because A^ = (ipi,Aipi), the left side of (21) is 
smaller than or equal to the right side. Therefore, only the opposite inequality is to be shown. Choose 
any orthonormal set <f>i, . . . ,<j> n of vectors of Ji and let (pi — J^JLi o-ij^j- Then 



n n oo n 

Ys^uMi) = EEi a «i 2A ; = E 


n oo 

+ E 


i—l i—1 j 


=1 i=l 


j=l J=n+1 


n 

* E 


n / n \ 

E m 2a ? + A «+i 1 - E i a «i 2 




_i=i 


V 3 = 1 / . 


n / n 

= nX n+1 +Y ( El^'l 
j=i \i=i 


(Aj - A„+i) < nA 



»=i j=i 



i=i 



which remains valid if we take the supremum in the leftmost member. 

Applying this lemma to A = a\(N) we see that X)/=o( n i) PH N < Sj=o A i f° r an y J- F rom (14), 
therefore, 

2 I™ A, (N)/N > 1 - - . (23) 

' * /V — >nn 



Lettig J tend to infinity we obtain (20). For Gross-Pitaevskii scaling in three dimensions Seiringer 
[17] proved the much stronger result limjv^oo X (N)/N = 1. 

The logic of the use of Lemma 2.1 for proving GBEC is somewhat different from the strategy 
followed in Theorem 1: First we make an appropriate choice of J and then we impose the necessary 
condition on the interaction. An obvious consequence of Lemma 2.1 is 

Proposition 2.1 Let J = J(N) with J = o(N) or J = [sN\. Choose 8 = £j+i — \{eq + E\). Fix 
(3 > and suppose that 



lim ■ 



N 



J 1 

E^'W&W) _ li( u N)fm% ~ mf U N ] 



3=0 



= b > 



(24) 



where lim means liniTv^oo if J/N — > and lim s ^ limAr^oo if J = \_sN\ . Then 

J 

A 



1 J 

lim AT E^'^ H « - 6 



(25) 



The inequality (25) implies at least GBEC. Note that 5 = S(J,N) with a separate dependence on N 
through the eigenvalues of H° if V is scaled. In some applications the positivity of the left member 
of (24) can hold with a single term, j = 0, of the sum. 

The results of the theorem and the proposition above apply to the ground state as well, if we take 
first the limit (3 — > oo, then N — > oo. In the use of Lemma 2.1 for proving BEC or GBEC the lower 
part of the spectrum of H° is shifted upwards in such a way that the ground state of the modified 
Hamiltonian H°(J,S) remains ipo, the ground state of H°. Therefore 



lim {n ) m o ( j S) = N 
and there is BEC or GBEC in the ground state of the interacting gas provided that 



(26) 



lim— [($ ,U N $ ) -inf U N ] < 1 . 
ivd 



(27) 
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It is interesting to see why and how the method described above fails in proving BEC or GBEC 
of the homogenous Bose gas. We exhibit this failure on the ground state. Let us consider first the 
one-dimensional homogenous 5-gas, that is, u(x) — 2c5(x) and 

V & = { +00 x < 0, x > L ^ 

or take a periodic boundary condition at and L. Let p — N/L. In the periodic case (<&o, Un&q) = 
cN(N — 1)/L and therefore GBEC follows from Proposition 2.1 and equation (27) if lim J/N = and 
cp < S ~ £j oc J 2 /L 2 = p 2 (J/N) 2 which means limc/p = \m\(J/N) 2 = 0, i.e. a vanishing interaction. 
Still, Proposition 2.1 yields a nontrivial estimate for interacting homogenous gases. Let, in general, 
Un — Sj<j u ( x i — x j) ^ —BN, J \u\ dx < 00, p > and consider TV bosons in a <i-dimensional 
cube of side L = (N/p) 1 ^. For periodic boundary conditions (<&o, Un&o) = \p{N — 1) Jwdx and 
therefore (27) holds true if 5 > \p J wdx + B > 0. Recalling that in the setup of the proposition 
lini/v^oo £j/S — 1 if J — > 00 with iV, choose 

J(JV) = [w d (2m<5) d / 2 /^ d ]iV (29) 

where Vd is the volume of the unit ball in d dimensions. Then the zero temperature version of 
Proposition 2.1 for homogenous gases proves 

1 J ^? ; 4 y fudx + s / 
J im i im n E ( n >)w» ^ 1 - iLJ —k > 30 

N^oo 8^00 iV ' — ' 

i.e. a macroscopic number of particles are distributed over a macroscopic number (29) of lowest lying 
levels. Although this does not prove BEC, the bound (30) is nontrivial because the number of levels 
is infinite for finite N. 

To use Proposition 2.1 at positive temperatures, BEC or GBEC has to be proven in the nonintcr- 
acting gas with the shifted spectrum. This is a minor problem if the external potential is fixed, as 
in the case of our forthcoming discussion of GBEC in superharmonic traps. For the scaled harmonic 
potential we will have to pay somewhat more attention to this question. 



2.2 Superharmonic traps in one dimension 

In one dimension the particularity of potentials with x 2 /V(x) — > as \x\ — > 00 is that n/e n — > as 
n — > 00. As an example, for homogenous potentials V(x) = c\x\ n semiclassical quantization yields the 
eigenvalues in the form 



htjc 



1/77 



in + 0(1)) 



Jo 



:dx 



r(f)r(iA?) 
r(| + i/v) ' 



(31) 



showing that the eigenvalues increase faster than linearly if 77 > 2. Note that Sturmian theory [27] 
confirms the semiclassical formula (31). 

Let, therefore, H° be a fixed one-particle Hamiltonian with a spectrum such that n/e n — ► (which 
is our definition of a superharmonic trap in one dimension) and choose J(N) — > 00 in such a way that 
lim J/N = and lim N/e.j = 0. Then for 5 = — ^(eo + £1) we also have lim N/ S = 0. 

In [14] we proved that in the case of a noninteracting gas in a fixed trap there exists a uniform 
upper bound on the mean value of N' = N — n , namely, for any N and any p < E\ — £q 



(AO 



< 



(1 _ e -^y 



ra=l 



1 

e -8(E„-£ -n) 



(32) 
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Similar inequality holds for (^V')/3i?^(J,5)i with the exception that the bound still depends on N because 
e n — eo has to be replaced by e n — Eq — <5 if n > J. However, due to superharmonicity, for n > J + 2 
we can use the estimate 

1 



En — £0 — $ — En — £,J+1 + "(El — £o) > E n -J-\ — Eq (33) 

if N (and thus J) is large enough. Choosing e.g. [i = \(e\ — £o) ; f° r sufficiently large N we obtain 



(N'}p H %(J,6) 



< 



(l-e-/»f)3 



n- 



1 



-i 2 



.ra=l 



e -0{e„-e o -lA 



Therefore 



l™^<no>/jflO (Ji5) 



= 1 



and the inequality (24) holds if 



(34) 
(35) 
(36) 



In words, GBEC follows if the mean interaction energy (Un)/3h° n does not exceed N 5 w Nej, the 
maximum energy of N noninteracting particles in the would-be (generalized) condensate. 
If Un is a pair interaction, we can use the estimate (12) and 



l(^,W.)|< ^ ^ IMIiIIpSIIi 



to obtain 



N -I, 



Polli + c (/^) Ihjvlli - 



inf C/at 
AT 



(37) 



(38) 



If the interaction is stable, inf Un/N > —B for some constant B; if it is not stable still inf Un/N > 
— (N — l)infuAf/2. In cither case, for un — u integrable, bounded from below and independent of 
N the quantity in the square bracket is of the order of TV and thus the limit in (38) vanishes. With 
Lemma 2.3 this implies 

Theorem 2 In one- dimensional superharmonic traps bosons interacting via a lower semibounded 
(unsealed) integrable pair interaction undergo a complete generalized Bose-Einstein condensation at 
all temperatures: For any choice of J(N) such that lim J/N — and lim N/ej = 



1 J 1 J 



(39) 



3=0 



The claim of the theorem could be made stronger. For instance, in a box (28) with a fixed L 
independent of N, J and S ~ Sj can grow almost as fast as TV and -/V 2 , respectively, and therefore we 
can obtain GBEC for interactions as strong as (Un) = o(N 3 ) instead of 0(N 2 ) stated in the theorem. 



3 Scaled harmonic trap in one dimension 

3.1 Bose condensation via phase transition in the noninteracting gas 



The one-particle Hamiltonian is 



„o V d 2 1 
2m dx 2 2 



2„2 



(40) 
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Measuring the energy from that of the ground state, the partition function for N noninteracting 
particles reads 

N 

Qnm = e- a ^m = e~ a Zm = J2q m ,a (41) 

{n J } j > :Y.ri J =N {nj } j>0 :£ ™j < N m =° 

with 

q m ,a= 53 e~ a Xw (42) 

{nj}j>o-T. rij=m 

where we have introduced the notation a = hw(3 and used f3(ej — So) = ja. The key to the forthcoming 
analysis is 

Lemma 3.1 

N 

QN,a=H(l-e- ka r 1 (43) 
fe=i 

m 

q^a^e-^Hil-e-^r 1 (44) 
k=l 

and the probability of having m particles in excited states is 

N 

PnA^' = m) = e- ma H (1 er ka ) . (45) 

k—ra+l 

Proof. 

PnA N ' =m) = q m ,a/QN,a and (44) follows from (43) through q m<a = Q m . a - Q m -i,a- Now Q N . a 
can be rewritten as 

oo oo oo 

Qn. a = 5Z e~ a ^=^ = ^ e- lia ^ e- t2a --- 53 e ~ lNa ( 46 ) 

0<ii<--<ijv ii— %2— ii iN—^N-i 

from which (43) follows. 

An alternative way is to compute first q m ,a by using that it is the generating function of p m (n), 
the number of (unordered) partitions of n into m parts, 

oc oo 

q m .a = 53 Pm(n)e- na = Y,Pm(n)e~ na (47) 

n—m n—1 

because p m (n) = for n < m. Starting with the identity 

Pm(n) =p m (n-m) +p m - 1 (n-m) H hpi(n-m) , (48) 

valid for n > m > 1, one can derive the recurrence relation 

1m,a = 1 _ e _ ma 53 «M *>.« = 1 (49) 

fe=0 

from which (44) follows by induction and (43) by (41) and (49). 
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Due to the simple form of Pn,o,(N' = m) we can obtain precise asymptotic results on the distri- 
bution of n in the case of different scalings. In what follows, we discuss the thermodynamics of the 
noninteracting gas in a scaled harmonic trap, characterized by the Hamiltonian 

N n 2 d 2 

H n = Et-^T^ +lNV{^ Xi )] (50) 



i=l 



where V(x) = ^mut 2 x 2 and jn — ► 0. This amounts to open the trap by replacing u> by uj-jn or, 
equivalently, a = fkj/3 by 

a n — a-fN — * 0. (51) 

From Lemma 3.1 



N 



P N ,a N (N' < m) = = I I" (1 - e~ la n . (52) 

l=m+l 

Let A^v = 1 — m/N where m is any integer between and N. Taking the logarithm of equation 
(52) and expanding the right member in Taylor series we find 



N 

lnP N , aN (^>X N )=-^A Nk , A NkM = - k £ ^- C ^Zrr iJ • ^ 

fe=l l=m+l 1 ' 

A-Nk(^N) are nonnegative, monotone decreasing with k and increasing with A^v- Later on, 

P N , aN Q > Xn) < e~ A "^ (54) 

will be used. 

The following two propositions serve to find necessary conditions on the sequence a» which give 
rise to a trivial asymptotic distribution of n /N. In this section we use the notation (no)N,a N for 

( n o)pH° N - 

Proposition 3.1 Suppose that 

X N -> A > and InJV- ^1 - Na N -> oo . (55) 

Then 

PN.a N ( ^7 > X N J -> and lim — (n )N,a N < X . (56) 

V iv / N^oo iv 

7/lniV — iVajy — > oo t/ien ^ (no)jv.ajv ~~ * 0; *- e - ^ere is no Bose-Einstein condensation. 
Proof. 

To obtain the vanishing probability it suffices to show that ^4jvi — ► oo. From the inequalities 

xe x/2 <e x -Kxe x (57) 

valid for x > 0, 

Awi > e lnW-(l-A w /2)Wo N +lnA w -ojv ^ ( 58 ) 

indeed. Moreover, 

-^■("-o)w,ajv < X N Pn,cl n (jfj: < Ajv) + PN.a N (J^r > X N ^ -> A. (59) 

If the stronger condition In TV — Na^ — > oo is fulfilled, (55) and therefore (56) hold true for all A > 0. 
This implies the absence of BEC. 

The counterpart of Proposition 3.1 is 



1 _ N _ p -kNa N I p k\ N Na N _ i \ 
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Proposition 3.2 Suppose that 

X N -> A < 1 and (1 - X N )Na N - InN -> oo . (60) 

T/ien _^ 

fW.ajv ( ^7 > A w ) -» 1 and Jim — (n ) N , aN > A . (61) 

V iV / iv — >oo iv 

If Ncln / \xiN — > oo i/ien -^{n a )N,a N — * 1; *- e - ^ere is a complete Bose-Einstein condensation. 
Proof. 

Using (57) and 

-ln(l-e-)<-^ T (x>0) 

we obtain 

oo 

< ^ Ajvfe < TVAat 
fe=i 

which implies the result on the limit of the probability. On the other hand, 

Jj{no)N,a N > AiV-PjV,a w > Ajv) -> A . (63) 

Suppose now that iVa at / In TV — > oo. One can choose a sequence A at — > 1 in such a way that (1 — 
X N )Na N /lnN -> oo. Then (60) is fulfilled and (63) holds for A = 1. 

From Propositions 3.1 and 3.2 we can identify the scalings which yield trivial limits: 

Corollary 3.1 If Nj N /\nN -> then T c = 0. If Nj N / In N -> oo t/ien T c = oo and t/iere is a 
complete BEC at all T < oo. 

Proof. 

In the first case Naiv/\nN — > and thus In N — Najy — > oo for any a < oo i.e. T > 0, and Proposition 
3.1 applies. In the second case Non/ In TV — > oo for any a > i.e. T < oo, and Proposition 3.2 applies. 

Thus, we have found that a phase transition can occur only if N-f N / In N has a finite nonvanishing 
limit. In this case without restricting generality we suppose that 

7A r = In N/N, a N = a In N/N. (64) 

Choosing a different prefactor can only change the critical temperature. Propositions 3.1 and 3.2 
already yield the phase transition: 

Corollary 3.2 If aN = alnN/N then a c = 1. For a < 1 there is no BEC. For a > 1 and < 

lim Xn < 1 — a -1 

j^ p — G^ A ") =L (65) 

Moreover, 

J im Tf( n o)w,a JV > 1 - - • (66) 
N^oo iv a 

Proo/. 

If a < 1, In AT — TVa^ = (1 — a) In AT — > oo and by Proposition 3.1 there is no BEC. If a > 1 and 
< A — lim An < 1 — a -1 then (60) and thus (61) hold true. Letting A tend to 1 — a -1 we find (66). 

This is a temporary result. We will show, among others, that for a > 1 the probability of having 
n /N > 1 — a -1 tends to zero and thus the distribution of uq/N becomes degenerate and concentrated 
onl-fl -1 . This will imply that in (66) there is equality. 



In 



-{l-\ N )Na N 



< 



A 



N 



>(1- \ N )Na N -In N _ 



(62) 
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Proposition 3.3 Suppose that 

\ N ->\<1 and (1 - \ N )Na N —> oo . (67) 

Then 

PN. AN (^>X N )=e-^ + ^ A ^^ (68) 
where < e N < - In (l - e -(T--*N)Na„j _^ . 
Proo/. 

Let us rewrite equation (53) in the form 



Now 



and therefore 



as N — * oo. 



lnP N , aN > A N ) = -A N1 [l + Y: • (69) 

Ajv.fc+1 e -(m+l)(fe+l)a N X e _;(fe+i )aN 

j4]V1 ~ e -(m+l)ajv fc + 1 e - (a « 

<- j; g-(m+l)fcaw ^ jj_ e -fe(l-Aj V )A r a N /•yg'v 

fc + 1 fc 

fe=l ^ Ari fc=l ^ 



Henceforth, we concentrate on the phase transition. With the scaling (64) the condition (67) is 
satisfied if a > and A < 1. Most of the results listed in the theorem below follow from Proposition 
3.3 applied to this case, 

/ran \ f ]ya\ N _ j ~| 

PN,a N > AjvJ = cxp |-iV-° ealnAf/jv _ i (l + e N )j. (72) 

Theorem 3 The scaling ajv = alnN/N leads to a phase transition at a = nw/3 = 1 no BPC for 
a < 1 and BEC for a > 1. Ira details, the following hold true. 
I. Limit distribution of hq . 

(i) For < a < 1 arad x > 

^ For a > 1 arad A^v — > A e]0, 1[ 

i- t-> / n o^> \ f 1 if A < 1 - a -1 ,_ 
n^oo aN \~N ) = { if A>l-a-. < 74 > 

fmj For a > 1 arad aray rea^ a; 

lim Pv. ajv Cr^T — 1 + a_1 l ln JV-o _1 IiiInJV > x) = cxp {-a-V x }. (75) 

N — >oo ' V L J\ J / 



Equivalently, 



^PN, aN (^(«o ~ («<>>*,„ J > = exp {-a-^^W) } 



(76) 
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where 



II. Mean value of no ■ 
(i) For < a < 1 



(ii) For a > 1 



/oo 
xexp {ax — a~ 1 e ax } dx 
-oo 



JV^oo N a 



TV 



(n ) N .a N = (1 - a^N + J^K 1 InlniV + r?(a)] + o 



TV 
InTV 



In particular, 



lim ^-(n )Ar, ajv = 1- -. 
at^oo TV a 



1/7. For any a > T/ie /ree energy of TV particles is 

N 6/3alnTV 



+ o 



TV 
In~TV 



(77) 

(78) 

(79) 
(80) 

(81) 



Proof. 

(Li) Substitute X N = x/N 1 ^ - into equation (72). Because A^lnTV — > 0, TV aAjv - 1 can be replaced 
by aAjvlnTV and (73) follows. 

(I.ii) If A < 1 — a -1 then Ani(\n) — > and if A > 1 — a -1 then Aat^Ajv) — > oo. Therefore the limit 
of (72) is the degenerate distribution (74). 

(I.iii) The Gumbcl distribution (75) is obtained by substituting 



1 In In TV + ax 
A at = 1 1 — 

a a in TV 



(82) 



into equation (72). Because A^y — ► 1 — a -1 < 1, condition (67) is satisfied and Proposition 3.3 indeed 
applies. The form (76) follows from (II. ii) to be shown below; 77(a) is the expectation value of the 
Gumbcl distribution. 
(Il.i) For any A > 

OO 

£ mA [p N , aN > mA) P N , aN ( J£ > (m + 1)a)^ 

< {J %^ < E ^ [P,, ON > (m - 1)A) - > ™A 



m—1 



(83) 



The sums are actually finite but the upper bounds tend to infinity with TV. Using the inequalities 
(54) and (57) 



P N .a N (^>x)< cxp{-A N1 (xN- 1+a )} < exp{-xe- alnN / N } < e^/ 2 



(84) 



for TV large enough. For the left and right members of (83) mAe"' m_1 ' A / 2 is a summable upper 
bound, thus we can interchange the limit TV — > 00 and the summation over m to find, with (73) and 
the convexity of the exponential, 



e- A J2 mAe- mA A < mA( e - mA - e -( m +^ A ) < lim 



m—1 



m—1 



JV— >oo 
oo 



("o)jV,aiy 
TV a 



< J2 mA{e- {m -^ A - e,- mA ) <e A J2 mAe~ mA A . 



(85) 



m—1 



m—1 
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Letting A tend to zero the Riemann sums go to J °° xe x Ax = 1. 
(II. ii) We start as before. Let 

f N (n )= + In AT -a" 1 In In TV. (86) 

For any A > 

oo 

mA [P N . aN (/jv(no) > mA) - P N . aN {f N {n ) > (m + 1)A)] 

m— — OO 

OO 

< (fN{n ))N, aN < mA t Pw ^« CM«o) > (m- 1)A) - P;v,a w (/jv(no) > mA)] . (87) 

m— — oo 

The sums are finite with the upper and lower bounds tending to infinity as N increases. Suppose that 
we can interchange the summation with the limit TV — > oo. Then (75) yields 

OO 

]T m A[e-»- le ™ A -e-»- le * + " 4 ] 

m— — oo 

oo 

< lim (f N (n )) N ,a N < V mA [e- a_lea< ^ 1)A - e-^ 1 ^! . (88) 

m— — oo 

Now cxp{— a~ 1 e ax } is concave if x < a In a and convex if x > a _1 lna. Accordingly, we divide 
the sums in two parts, bound the differences in the square brackets with A times the derivatives 
at the upper or lower end of the intervals and let A go to zero. Both the upper and lower bound 
°f (/iv(^o))iv,aiv converge to 77(a). Equation (79) follows by simple rearrangement. In (87) the in- 
terchange of the summation with the limit N — > 00 is again based on the dominated convergence 
theorem. However, the sums have to be divided in two parts. For m > we can use (54) with 
(82) and x = (m — 1)A as an upper bound on the difference (actually on both terms) in the square 
brackets because Ajvi(Ajv) > (2a)~ 1 e ax if N is large enough. For m < only the difference in the 
square bracket is small. Using (71) a lengthy but straightforward computation yields 

P N , aN (/jv(no) > mA) - P N , aN (f N (n ) > (m + 1)A) < 1 - cxp {-2 a -\e aA - l)e amA } (89) 

which decays exponentially as m — > — 00 thus yielding a summable upper bound. 

Equation (80) is a consequence of (79) but can also be obtained directly from the degenerate 
distribution (74): Choose any e > 0. 



(1- a- 1 - e) P N , aN Q >1- a" 1 - e) < 



(no) 



N,a N 



N 



< (1 - a- 1 + e) P N>aN < 1 - a" 1 + e) + P N . aN > 1 - a" 1 + e) . (90) 
Taking the limit N — > 00 and applying (74), 

1 — a — e<hmmf ■ — < hmsup ■ — <1 — a +e 91 

which holds for any e > and thus for e = as well. 

III. If Zx, aN denotes the TV-particle partition function then 

F% = -^In Z N , aN =~ (-^Na N +lnQ N , aN ^ . (92) 
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From (43) 



N 



mQ N , ajv =-]Tln(l- e - fcQ «) (93) 
fe=i 

so that 



lim 7iY lnQw.a^ =- lim V 7iY In(l - e~ ak ^) = - / ln(l - e~ ax ) Ax (94) 

N^oo N^oo * — J Jq 

because jn — * and A^jv = InAf — ► oo. Expanding the logarithm and integrating term by term 
Efe=i e~ kax /k — ► — ln(l — e~ ax ) monotonically) we find n 2 /6a and, hence, equation (81). 

(1) Equation (76) implies |no — {no)N,a N \ ~ AT/lnAf for a > 1. These are huge fluctuations even 
compared with the super- normal fluctuations |no — (no)| ~ A^ 2 / 3 in the condensation regime of the 
three-dimensional homogenous Bose gas [28] . 

(2) There is no singularity in lim Fn = ~ Tr 2 /6/3a at the critical point a = 1. 

(3) For 7jv = 1/-/V equation (94) yields an extensive free energy 

lim = -/T 1 lim 1 \nQ N , aN - / ln(l - e~ ax ) dx (95) 

so that this scaling corresponds to the homogenous limit. According to Corollary 3.1, T c = in this 
case as it has to be in the one dimensional homogenous Bose gas. 



3.2 Generalized Bose condensation at all temperatures in the noninteract- 
ing gas 

In the three-dimensional homogenous noninteracting Bose gas above the critical temperature the 
mean occupation number of each one-particle state remains finite in the thermodynamic limit. This 
is immediately seen in the grand-canonical ensemble and is valid in the canonical ensemble as well, 
due to the strong equivalence of ensembles (the Kac density is a Dirac delta) . The situation is quite 
different in the one-dimensional scaled harmonic trap when = I11N/N. As we shall see, there is 
a complete GBEC at all temperatures. It will also be shown that below the critical temperature the 
condensate is not fragmented, (ni) N ^ aN — o(N) for i > 0. This means that for a > 1 a condensate in 
the ground state of H° whose density is 1 — 1/a coexists with a generalized condensate of density 1/a. 

The intuition behind the results of this section is guided by the following observation. (We drop 
the subscript on which plays no role here.) 

Lemma 3.2 For any i 

N / i-1 \ 

5>,- = M . (96) 

M=0 \j=0 J 

Proof. The reader can easily check that because of 

e n = e + n{ei - e ) 

the conditional distribution of n i; given the number of particles in the lower lying eigenstates, satisfies 



Pn rii 



i-1 

Y^rij =M I =P N _ M {n =m) . (97) 
3=0 
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Multiplying by m and summing over it the right member becomes (no)iv-M while on the left-hand 
side we obtain the conditional expectation value (rii\no + • • • + rij_i = M)jy. Multiplying by the 
probability of the condition and summing over M yields (96). 

Because for a < 1 n ~ N a , we expect that 

(n t )N,a N ~ (n )N-iN°,a N ~ N a (98) 

if i <C A^ 1_a , and that an o(N) number of lowest lying levels carry roughly all the particles. This 
should hold true for a > 1 as well and, because of n ~ N(l — 1/a), we expect also that 

{ni)N.a N ~ (no)iv/a-(i-i)a- 1 ivininJV/iniv,o N ~ a ~ 1 N In In AT/ In N (99) 

for 1 < i <C In A^/ In In N. Not all these conjectures will be verified below. We start by proving the 
second part of (98). 

Proposition 3.4 All the results of Theorem 3 remain valid if the scaling ajy = a In N/N is replaced 
by a' N = (1 + t}n)o,n where ijn = o(l/ In N). 

Proof. Consider A^ki defined in (53), as a function of on and Ajy. If a/v is replaced by a' N given 
above, one finds that for any sequence A at 

ANk(a-N, ^N)/Aiyi c (a' N , \n) — * 1 
at least as fast as tjn^N tends to zero. Therefore equation (68) can be replaced by 

P N ,a' N (^ > A W ) = e -(l+C,)A N1 (a N ,X N ) (10Q) 

where — > 0. Because the results of Theorem 3 grouped under points I and II were derived from 
equation (68) and did not depend on the particular way cm tended to zero, the modified scaling will 
provide the same outcome. Although we do not use it, we note that the free energy also will be the 
same. 

As a corollary we obtain 
Corollary 3.3 Let S = S{N) = o(N/\nN). For a < 1 

Urn N- a (n ) N +s,a N = 1. (101) 

JV^oo 

Proof. Because ]xm(N + S) a /N a = 1, 

lim N- a (n ) N+s ,a N = lim N- a {n ) N ,a> 

with 

1 + (lnA^)- 1 ln(l - S/N) . . . . 

° !n = 1- S/N a N = (l + VN )a N . (102) 

One can easily verify that tjnIiiN — > 0. Thus, the result follows by applying Proposition 3.4. 

Lemma 3.3 For any noninteracting Bose gas with a one-particle spectrum {sj} the average occupa- 
tion numbers in the N -particle canonical ensemble are strictly decreasing with the energy, 

(rij} N < (ni) N if Si < ej. (103) 
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Proof. Let Xi = cxp(— /3si) and let Zm and Z^^j denote the TV-particle partition functions of the full 
system and of the system with missing levels i and j, respectively. 



Z N [{ni) N - (nj) N } = J2( m ~ l)x?x l j Z N-i- m ,i,j = ( TO ~ l )i x T x j ~ 4 x ?] Z N-l-m,i,3 (104) 

l,m 0<Km<N 

= (m- l)(x iXj ) l [x?- 1 - x?- l ]Z N + miitj > 0. (105) 



0<Km<N 



Theorem 4 For any a = fruifi > the scaling a^v = a\nN/N leads to a complete generalized Bose- 
Einstein condensation. For a > a c = 1 there is no fragmented condensation but the generalized 
condensate on the levels 1,2,... is at its critical point. In particular: 
(i) For any a > 



(n) For a < 1 

(m) For a > 1 
but 



j<2/a N 



Jim AT a (^)jv,a« = 1 if i = o 



N— >oo 



iV 



1-a 



In AT 



lim AT ^n^Ar.ajv =0 if i > 1 



lim AT ''(n^jv.aw = 00 / or an 2/ ?? < 1. 

AT— >oo 



Proo/. 

(i) Let J = 2/qn- We shall prove that 



In P; 



N,a N 



=rnj < - (2- m if 



m> J 



Then for & = e - 2+7r2 / 12 



Y mP N,a N l^nj = m\ < Y m&m < 

m>J \j>J J J<m<N 

as N tends to infinity. Thus for A^ large enough 



J 



1-6 (1-6) 2 



(106) 

(107) 

(108) 
(109) 



(110) 



(111) 



0<N- ^2{nj) N , aN = ^(n^jv.ajv = mP N,a N 



n,- = to 



kj>j 



+ Y mP N,a N \Y n '<= m ) < J + 1 ' 

m>J \j>/ / 

Dividing by N and taking the limit we obtain (106). 

To prove (110), consider the probability on the left side. With the notation x = e~ aN 



<3>J 



\{ n j}j>J : T, n } =m 



( 



(112) 



• (H3) 
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The quantity in the first bracket is x mJ Q m .a N , cf. equation (41). The quantity in the second bracket 
can be bounded by dropping the constraint n j < TV — to, 

e ^ ,jnj I < n t^i = ( iw ) 

\{« 3 }/=M: n rf <JV-m 

see equation (43). Thus, 



PN.a N \^ n 3 = mj < x mJ Q N ) aN Q m , aN Qj- haN < x mJ Q, ; ^ aN (115) 
because Q m ,a N is increasing with m. Taking the logarithm and using Jcln = 2, 

\nP N , aN Ij2 n 3= m ) < -2m- ^]ln(l -x J ). (116) 
Now — ln(l — x- 7 ) is a positive decreasing function of j, therefore 

J-l , ,2 

- Vln(l-e- jaN ) < / In(l-e-*)dt 

p[ a N Jo 

< -- / ln(l - e"*)dt = Jtt 2 /12 < tott 2 /12 (117) 
2 Jo 

for m > J. Sustituting this into equation (116) we obtain (110). 
(ii) To prove (107) we remark that because of 

(ni}N.a N < (no)N,a N ~ N a , 

cf. Lemma 3.3, it suffices to show that lim /V^^at^ > 1. Let K = i In TV; then K = o{N 1 - a / In TV). 
From (96) 

/ i-i \ / i-i 

(n t ) N , aN > ^2 ( n o)N-M,a N PN,a N I E ^ = M - ./V- L-fS"JV a J ,a N P N,a N ^«j< ifTV a 

M=0 \j=0 / \j=0 

(118) 

where we have used the monotonic increase of (no)k,a N with fc, derived in [29]. By Corollary 3.3 
equation (118) implies 



lim iV-^iv,™ > 1 - lim P N>a „ W-> ^/v a . (119) 

N— >oo iV— >-oo \ z — ' ' 



We show that the limit on the right-hand side is zero. 



p N,a N \ Yl n J ^ L I = P N.a N ( -Ynj > - ) < p N.a N (max{n 1 } > -J 

< E (»i >y)=E E * roj E - E fe " fc 

J=0 V / j=0 m>L/i {n k } k ^j-.Y^ n k =N-m 
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i-1 i-1 

< xL3 ' tp N.a N (n > L/i, n 3 = 0) < PjV,a w ("o > S xLj/i 

3=0 3=0 



For L = KN a = iN a In TV we can apply (72) with X N = InN/N 1 - 11 . We find 

^ e aln 2 A7AT~ 



(120) 



Pn,o, n 

\j2 n J> KNa ] < ex p 



]\fa ga In N/N 



.Lj/i 



3=0 



(121) 



which proves the assertion. 

(iii) Choose any S such that 1 < S < a. Let if = (S/a)N and for the sake of notational simplicity 
suppose that if is integer. From equation (96) 



< (no)K,a N PN,a N (no > N - if) + NP Ny a N (n < N - if). 

Now N = (a/S)K, thus 

lnif + ln(a/<5) , nx ,MK 
if N is large enough. Because (no)K, a is an increasing function of a (or of (3, see [29]), 



(122) 



1 6 

< lim — (ni)jv iOW < - 
N— >oo iv a 



1 



1 



2(5-1 



lim Pjv.ajv (n < N - if) 

N— >00 



1 



1 



2(5 - 1 



+ 1 - lim P N „ N — > 1 - - 



N^oo 



n 



N 



a J a 



1 



2(5- 1 



(123) 



where we have used equations (80) and (74). Because this holds for any S > 1, letting S tend to 1 we 
obtain (108) for i = 1. For i > 1 (108) results by applying Lemma 3.3. 
To prove (109) choose any r\ < 1. From equation (96) 

(ni) N .a N > ( n o)M,a N PN,a N (no = N - M) > (n ) K , arr [ 1 - p N,a N (n > N - if)]. (124) 

M>K 

Let 5 = |(T7 + 1), then <5 < 1. Define X = (5/o)JV. By (74) now 

lim P N , aN (n > N-K) = lim P N , aN (n /N > 1 - 5/a) = 0. 

N^oo TV^oo 



Because > <5 In if/ if, with equation (78) we obtain 



Jim T7r("i)w,a N > 



sy 

- lim if 

a / AT— >oo 



1-5 



1 

if* 



5 In K/K 



X\ r > 

lim if 1 " 5 = oo (125) 
a / K^oo 



which finishes the proof of the theorem. We note that with more effort, using equations (76) and (79) 
one could prove the precise asymptotics (99) for, at least, any fixed i > 1. 



3.3 Generalized Bose condensation at all temperatures in the interacting 
gas 

In this section we investigate the possibility of introducing a nontrivial interaction in the system 
without loosing Bose-Einstcin condensation. By nontrivial we mean an interaction providing an energy 
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contribution comparable with or larger than the free energy of the free system. The Hamiltonian of 
this latter is (50) with — \nN/N while that of the interacting gas is 



Hn — H% + Un — ^ 



N r fi 2 d 2 



i=l 



+ U N (x u ...,x N ). (126) 



We return to the notations introduced in Section 2.1. Mean values with respect to the noninteracting 
gas in the scaled harmonic trap will thus be denoted by (^pn^i in contrast to (-}N,a N used in Sections 
3.1 and 3.2. Also, in the case of the spectrally deformed noninteracting gas we apply the notation 

(■)()H° N (J,6)- 

Theorem 5 Interacting bosons in a scaled harmonic trap with 

7Ar = and U N > —BN, {U n ) 0h o n = o(N\nN) 

undergo a complete generalized Bose-Einstein condensation at all temperatures: For any (3 > 

sN 

i-J™ iv&^ = 1 - (127) 

3=0 

For stable pair interactions, 

U N {xi, . . -,x N ) = y^u N (xj - Xj), 

i<j 

{Un)[)h q n — o(N\nN) and thus (127) holds true if ||mjv||i = °{in)- If u n( x ) — b N u{a^x), the 
condition reads \\u\\\ < oo and bx = o(a at In iV/iV). 

Observe that the scaling condition ||mat||i = o(\nN/N) is somewhat weaker than (3) and, again, 
examples from mean-field type to sharply concentrated interactions can be obtained. 

We note that the discussion of Section 2.2 applies to the case of a fixed harmonic potential. In 
analogy with Theorem 2 we find that for a fixed harmonic trap and scaled interactions satisfying 
{Un)ph°, = o(N 2 ) there is a complete GBEC at all (3 > 0. In Theorem 5 we conclude that for 
somewhat weaker interactions the same holds true even if the frequency of the confining harmonic 
potential is scaled with 7jv = hiN/N. One could also show without much further ado that for more 
general scalings lo n = 7ato;, where 7^ > In N/N, and (Un)ph% = °(N 2, Jn) there is a complete GBEC 
at all temperatures. 

Proof. 

We use Lemma 2.1 and Proposition 2.1. If J = |_ S ^J then 5 = e,j = fkuj]y([sN\ + |) rj sfkulnN 
implying 

lim m ^ Un ^ h °n _ inf Un ^ = 

and therefore 

1 3 1 J 

linl n X^'W" > I™ ^2(n 3 )f 3 H%(j,S) ■ 

3=0 3=0 

We can conclude by showing that 

( n j)l3H° N (J,S) > ( n j)^H% (128) 

because from Theorem 4 it follows that 

1 J 

lim — y^(rij}i gH a =1 

AT^OO _/V ^ 2 ^ N 

3=0 
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for any 0. To prove (128) notice that the spectrum of H°(J, S) — Eq — e,j is a part of the spectrum of 
\{H° -eo). Indeed, 

(129) 



spec ji(#°- £o )J = ± lN hwx {0,1,2,...} 



while 



spec {H°(J, 5) -e - ej} = ]p N two x {0, 1, 2, ... , 2 J, 2 J + 1, 2J + 3, 2J + 5, . . .}. (130) 



In [29] we proved the following result. 

Lemma 3.4 Let H° and H 1 be two one-particle Hamiltonians, both e~@ H and e~P H trace class and 

spec H° C spec H 1 

where repeated eigenvalues are considered separately. If is a common eigenvalue then the averages 
of the occupation number rn in the two noninteracting ensembles satisfy the inequality 

KWj, < KW° • ( 131 ) 

On the left- and right-hand sides of (131) denotes n[<pj] and n[<£>°], respectively, where H°tp9 — 
i?Vi = Ei<p}, and ip\ and may be different. If rjj are the eigenvalues of H 1 not contained in the 
spectrum of H° and Xj — e~^ Vi then 

N N 

teWj, - w^rr E K)^_ m ^[/3<- m ] E IK 3 = E K)^r m p ™ ( 132 ) 

^ Ni m=0 {fy} : E fe j= m J m=0 

where Z denotes partition functions. Because X)m=o ^™ = 1> (131) follows from the strict monotonic 
increase of {n^p^ with N, proven in [29]. We obtain (128) and thereby the proof of a complete 

GBEC by replacing H° with H°(J,5) - e - £j and H 1 with \{H° - e ) in Lemma 3.4. (To avoid 
confusion we precise that on the right-hand side of equation (132) H^_ m is defined with the same 
one-particle spectrum for each m.) 
In the case of pair interactions 

2 

{un) ph % = N , N _ ^ {Un)ph% (133) 
and to satisfy the condition of the theorem we need (ujv)/3ijo = o(lnN/N). We show below that 

K«JV>/jflo|<^(l + 0(V7KF))||«Jv||i (134) 

A/3 

and therefore ||ujv||i = °(jn) is a sufficient condition for GBEC to hold for any > 0. Now 

N 



(u n }(3h° n = Z\J3H%] 1 ^ / dx i dx 2 u N (xi - x 2 ) I dx 3 ■ ■ . dxN J\(xj\e 0H °\x g (j)) 

geS N J J j=i 

„ Tr £/(g)e _/3ff « f 

E ^ T T77r\ — bSS- / dx 1 dx 2 u N (x 1 - x 2 )ii g , l3 (x 1 ,x 2 ) (135) 



ge.S N 



where Sn is the group of the permutations of 1, 2, . . . , N, U(g) is the unitary representation of the 
permutation g in the full iV-particle Hilbert space and Tr is the trace in this space. In the average 
above there are two kinds of probability measures /U s ,/3- If 1 and 2 arc in different cycles of g then 

» 9 A X > V) = WiW^^l^ls)^!/) = M/3i {x)»fa (v) (136) 
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where £i(g) is the length of the cycle of g containing i and 

. . (x\e-P H °\x) 

If 1 and 2 are in the same cycle of length I and g^(\) = 2 then 

- jf3H °\y)(y\e-V-M H °\x) 



fj- g ,p(x,y) = 



{x e 



(137) 



(138) 



tle -t0HO 

In the first case by Fourier transforming and applying Schwarz inequality and Parseval formula we 
get 

J dx 1 dx 2 u N (x 1 - x 2 )fJ, g ,0{x 1 ,x 2 ) < ||mw||i||M/3iII ||m#J ( 13 9) 
with \\fip\\ denoting the usual L 2 norm of fj,p. In the second case 

J dx 1 dx 2 u N (x 1 - x 2 )fig,p(x 1 ,x 2 ) < J dz\u N (z)\ J dy fi g ,p(z + y 7 y) 

< \\unWxSup J dyfj, g ,p(z + y,y). (140) 

At this point we recall that (6) involves a transformation between two Hilbert spaces, that of L 2 
functions of the variables Xj and y, = ajyXi, respectively. Previously we have only been interested in 
properties depending on the spectrum of j3H°. Because (3 and w appeared in a single dimcnsionless 
combination a — hu/3, jn could be considered to multiply a and yield = ajN. Here we are 
estimating functions and have to decide which space to work in. In accordance with equations (50) 
and (126) we choose to scale the potential which amounts to replace u by uom — <^7jv and a by 
a at = and to keep (3 as an independent unsealed variable. 

We use Mehler formula [30] 



(x\e-^°\y) 



rriLUN 



to obtain 
Because 



27rfi,sinh ajy 

Hp{x) — \J s/7TCXp( — SX 2 ) 

dnp(x 



exp 



itiujn{x 2 + y 2 ) mujNxy 



2fttanhaAr /isinhajy 



s = s{(3) = — - — tanhftWArp/2 . 



ds 



Hp(x)[(2s) 1 - x 2 ] = jip{x)[{x 2 ) w - x 2 } , 



(141) 

(142) 
(143) 



\np\\ is an increasing function of s and, therefore, of (3: 



d_ 

ds 



[ l$dx = 2/ fX(3 ( x )?£l( x )dx + 2 [ ^ x )^l( x )dx 

J Jx 2 <l/ 2 s OS Jx 2 >l/ 2 s OS 

7 ^(*)cl*+ / ^dx 



2 <i/: 
> 2^{^/Tj2s) 



Here we used that np(x) is a decreasing function of \x\ and J /jpdx = 1. Thus, we can bound (139) 
by inserting the largest possible (3\ and (3 2 on the right-hand side. Because 



max£i( ff )/3 = (N - l)(3 < N/3 , 

9 
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wc find 



= [ e-^ N ^ 2 Ax = v^p = K tanh tuv N N/3/2 

n J L 2nn 



< 



mujN 



2irh 



1/2 



To estimate the right-hand side of (140), observe that 

v{z) = J n 9t0 {z + y,y)dy 



1/2 



(144) 



is normalized. One then obtains 



U(Z) = y/D/ 



ire 



D = 



mujN 



1 



+ 



1 



+ — 



1 



+ — 



1 



Ah \tanhjajv tanh(£ — j)a,N sinhjajv sinh(^ — j)a,N 



and thus sup 2 v(z) = u(0) = y/D/n. Now D attains its maximum for t = 2, 



(145) 



max I? 

e>2 



+ 



1 



2h \ tanh a» sinha^v 
Thus for the permutations containing 1 and 2 in the same cycle we have the uniform upper bound 

V2 



/v 2 



(146) 



Finally, wc obtain the inequality (134) by substituting (144) and (146) into (139) and (140), respec- 
tively, and using these latter to bound the integrals in (135). The consequence 6at = o(aN"fN) for 
scaled pair interactions is trivial. By this we finished the proof of the theorem. 



4 Summary 

This paper is an extension of our earlier study [14] of Bose-Einstein condensation of trapped Bose 
gases. We have concentrated on two problems. The first, treated in Section 2, concerns bosons 
in a fixed trap interacting with unsealed interactions; that is, neither the confining potential nor 
the pair interaction between two particles depends on the number of particles. Harmonic or weaker 
trap potentials combined with unsealed interactions are too difficult to deal with. However, if the 
confining potential increases faster than quadratically at infinity, in one dimension wc can prove a 
complete generalized Bose-Einstein condensation at all temperatures. This means that as N tends 
to infinity all but a vanishing fraction of particles will be distributed over a set of one-particle states 
whose number is asymptotically negligible compared with N. The result does not provide a more 
precise information about the number of one-particle states carrying the condensate. The generalized 
condensation may eventually prove to be normal and non-fragmented. 

The second problem we have investigated in this paper is the condensation of bosons in a one- 
dimensional scaled harmonic trap. In a previous work [25] it was shown that for a particular scaling, 
when the oscillator frequency u is replaced by uu In N/N, Bose-Einstein condensation occurs with a 
phase transition in the noninteracting gas. Sections 3.1 and 3.2 have been devoted to a detailed study 
of this system. Our most interesting finding is that the occupation of the excited states is anomalously 
large and results in a complete generalized Bose-Einstein condensation at all temperatures, superim- 
posed on the normal Bose-condensation below the critical temperature. In Section 3.3 we have studied 
the same system in the case when there is also a suitably scaled interaction among the particles. We 
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have shown that the complete GBEC is preserved by the interaction without being able to prove that 
the phase transition also persists. 

Acknowledgment 

This work was partially supported by the Hungarian Scientific Research Fund (OTKA) Grants T42914 
and T46129 and the Center of Excellence Grant ICA l-CT-2000-70029. 



References 

[1] M. Girardeau: Relationship between systems of impenetrable bosons and fermions in one dimen- 
sion. J. Math. Phys. 1 516-523 (1960) 

[2] T. D. Schultz: Note on the one- dimensional gas of impenetrable point-particle bosons. J. Math. 
Phys. 4 666-671 (1963) 

[3] A. Lenard: Momentum distribution in the ground state of the one- dimensional system of impen- 
etrable bosons. J. Math. Phys. 5 624-637 (1964) 

[4] P. C. Hohenberg: Existence of long-range order in one and two dimensions. Phys. Rev. 158 
383-386 (1967) 

[5] M. Bouziane and Ph. Martin: Bogoliubov inequality for unbounded operators and the Bose gas. 
J. Math. Phys. 17 1848-1851 (1976) 

[6] M. D. Girardeau: Broken symmetry and generalized Bose condensation in restricted geometries. 
J. Math. Phys. 10 993-998 (1969) 

[7] N. N. Bogoliubov: Quasi- averages in problems of statistical mechanics. Dubna report D-781 
(1961), in Russian 

[8] E. H. Licb and W. Liniger: Exact analysis of an interacting Bose gas. I. The general solution of 
the ground state. Phys. Rev. 130 1605-1616 (1963) 

[9] F. D. M. Haldane: Effective harmonic-fluid approach to low-energy properties of one-dimensional 
quantum fluids. Phys. Rev. Lett. 47 1840-1843 (1981) 

[10] D. B. Creamer, H. B. Thacker and D. Wilkinson: A study of correlation functions for the delta- 
function Bose gas. Physica 20D 155-186 (1986) 

[11] V. E. Korcpin, N. M. Bogoliubov and A. G. Izergin: Quantum inverse scattering method and 
correlation functions. Cambridge University Press (1993) Ch. XVIII. 2 

[12] S. Coleman: There are no Goldstone bosons in two dimensions. Commun. Math. Phys. 31 259-264 
(1973) 

[13] L. Pitaevskii and S. Stringari: Uncertainty principle, quantum fluctuations, and broken symme- 
tries. J. Low Temp. Phys. 85 377-388 (1991) 

[14] A. Siito: Thermodynamic limit and proof of condensation for trapped bosons. J. Stat. Phys. 112 
375-396 (2003) 

[15] K. Symanzik: Proof and refinements of an inequality of Feynman. J. Math. Phys. 6 1155-1156 
(1965), 

[16] E. H. Licb and R. Sciringcr: Proof of Bose-Einstein condensation for dilute trapped gases. Phys. 
Rev. Lett. 88 170409 (2002) 



25 



[17] R. Seiringer: Ground state asymptotics of a dilute, rotating gas. J. Phys. A: Math. Gen. 36 
9755-9778 (2003) 

[18] F. Dalfovo, S. Giorgini, L. P. Pitaevskii and S. Stringari: Theory of Bose-Einstein condensation 
in trapped gases. Rev. Mod. Phys. 71 463-512 (1999) 

[19] E. H. Lieb, R. Seiringer and J. Yngvason: Bosons in a trap: A rigorous derivation of the Gross- 
Pitaevskii energy functional. Phys. Rev. A 61 043602 (2000) 

[20] A. Gorlitz, J. M. Vogels, A. E. Leanhardt, C. Raman, T. L. Gustavson, J. R. Abo-Shaeer, A. 
P. Chikkatur, S. Gupta, S. Inouye, T. Rosenband and W. Ketterle: Realization of Bose-Einstein 
condensates in lower dimensions. Phys. Rev. Lett. 87 130402 (2001) 

[21] P. Bouyer, J. H. Thywissen, F. Gerbier, M. Hugbart, S. Richard, J. Retter and A. Aspect: One- 
dimensional behaviour of elongated Bose-Einstein condensates. Quantum gases in low dimensions. 
Les Houches Proceedings (2003) 

[22] E. H. Lieb, R. Seiringer and J. Yngvason: One-dimensional bosons in three-dimensional traps. 
Phys. Rev. Lett. 91 150401 (2003) 

[23] D. S. Petrov, G. V. Shlyapnikov and J. T. M. Walraven: Regimes of quantum degeneracy in 
trapped ID gases. Phys. Rev. Lett. 85 3745-3749 (2000) 

[24] F. Gerbier, J. H. Thywissen, S. Richard, M. Hugbart, P. Bouyer and A. Aspect: Momentum 
distribution and correlation function of quasicondensates in elongated traps. Phys. Rev. A 67 
051602(R) (2003) 

[25] W. Ketterle and N. J. van Druten: Bose-Einstein condensation of a finite number of particles 
trapped in one or three dimensions. Phys. Rev. A 54 656-660 (1996) 

[26] V. A. Zagrebnov and J.-B. Bru: The Bogoliubov model of weakly imperfect Bose gases. Phys. 
Rep. 350 291-434 (2001) 

[27] E. C. Titchmarsch: Eigenfunction expansions I. Oxford at the Clarendon Press (1962) Ch. VII 

[28] E. Buffet and J. Pule: Fluctuation properties of the imperfect Bose gas. J. Math. Phys. 24 1608- 
1616 (1983) 

[29] A. Siito: Correlation inequalities for noninteracting Bose gases. J. Phys. A: Math. Gen. 37 615- 
621 (2004) 

[30] G. Rocpstorff: Path Integral Approach to Quantum Physics (Springer, Berlin, 1994) 



26 



